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INTRODUCTION 


G.  Weller 

Weather  and  climate  in  the  arctic  and  sub-arctic  regions  are  strongly 
affected  by  the  state  of  the  Arctic  Ocean  and  the  atmosphere,  and  by  the 
effective  coupling  between  them.  The  polar  regions  occupy  the  role  of  the 
heat  sink,  as  part  of  the  global  heat  engine  which  drives  the  general  cir¬ 
culation  -  studies  of  these  heat  sinks  ar~i  not  only  important  in  under¬ 
standing  weather  and  climate  on  a  regional  scale,  but  in  fact,  also  on 
a  global  scale. 

In  the  Arctic  Basin,  two  important  components  of  the  environment 
strongly  influence  the  effectiveness  of  the  heat  sink,  which  owes  its 
existence  to  net  radiative  energy  losses  -  these  are  the  presence  of 
perennial  pack  ice  and  low  stratus  cloud  decks  during  summer.  ARPA- 
sponsored  research  by  the  University  of  Alaska  into  the  radiation-related 
properties  of  these  two  media,  and  into  the  general  solar-terrestrial 
radiation  matrix  in  the  Arctic,  commenced  during  spring  1971.  These 
studies  aim  to  aid  in  the  understanding  of  the  radiation  phenomenon  in 
the  Arctic  atmosphere  and  at  the  terrestrial  Arctic  surface,  and  ultimately 
relate  to  considerations  of  the  heat  budget  of  the  polar  regions  and  to 
climate  modifications  and  controls  in  that  area.  Our  studies  to  date 
have  concentrated  on  determining  the  physical  characteristics  of  both 
pack  ice  and  clouds,  as  they  affect  the  Arctic  radiation  regime.  At 
Barrow,  flights  through  stratus  clouds  with  aircraft-mounted  continuous 
cloud  particle  samplers  have  determined  number  and  mass  densities,  particle 
size  distributions  and  ice  crystal  contents  of  these  clouds  over  a  wide 
temperature  range.  Simultaneously,  light  intensity  and  reflectivity 
measurements  have  been  carried  out  to  determine  the  optical  properties 
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of  these  clouds.  A  study  of  precipitation  mechanisms  is  also  underway. 

Aircraft-mounted  photometers  have  provided  measurements  that  will 
allow  deduction  of  aerosol  profiles  up  to  13,000  feet;  on  two  occasions, 
these  were  extended  up  to  35,000  feet  aboard  flights  by  the  NASA  CV-990 

i 

aircraft  to  the  Arctic  Ice  Dynamics  Joint  Experiment  (AIDJEX)  camp  in 
i  the  Beaufort  Sea.  Although  the  air  is  fairly  "clean"  over  the  Arctic 

Ocean,  ice  crystal  aerosols  are  present,  and  are  important  in  scattering 
incoming  radiation.  In  addition  to  measurements  of  ice  nuclei  concentra¬ 
tion  at  Barrow,  actinometric  observations  through  frequently  thick  ice 
crystal  clouds  advected  from  open  leads  were  made  on  numerous  occasions 
on  the  AIDJEX  ciiinp. 

Routine  monitoring  of  components  of  the  radiation  balance  have  been 
carried  out  at  Barrow,  Ice  Island  T-3  and  the  AIDJEX.  camp.  Coupled  with 
%  this  were  studies  of  surface  inhomogeneities,  such  as  hummocks,  snow 

dunes,  leads  and  pressure  ridges,  and  the  effects  of  these  features  on 
the  surface  albedo  and  temperature.  Of  particular  interest  were  the  mea¬ 
surements  over  refreezing  leads,  both  artificially-made  and  natural  ones, 
showing  extraordinary  modification  of  the  radiation  balance  during  the 
freezing  process.  Light  extinction  measurements  in  snow,  ice  and  sea 
water,  using  photocells,  allowed  the  construction  of  a  typical  extinction 
nomogram  for  these  substances. 

To  synthesize  all  observations,  theoretical  and  computer  modeling 
is  playing  an  increasingly  important  part.  We  hope  to  pursue  this  latter 
'  aspect  much  more  vigorously  in  the  future. 

The  authors  of  this  report  are: 

Dr.  Sue  Ann  Bowling,  Assistant  Professor  of  Geophysics 
Dr.  K.O.L.F.  Jayaweera,  Assistant  Professor  of  Geophysics 
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Dr.  Takeshi  Ohtake,  Associate  Professor  of  Geophysics 
Mr.  Stan  Parker,  Research  Assistant 
Dr.  Glenn  Shaw,  Assistant  Professor  of  Geophysics 
Dr.  Gunter  Weller,  Associate  Professor  of  Geophysics 
Dr.  Gerd  Wendler,  Associate  Professor  of  Geophysics 

All  are  staff  members  of  the  Geophysical  Institute  of  the  University  of 
Alaska . 

The  results  presented  in  this  report  are  incomplete.  The  final  report 
on  this  contract  will  be  available  by  spring  1973. 


1.  RADIATION  REGIME  OF  THE  ARCTIC  BASIN  -  STUDIES 
AT  T-3,  AIDJEX  SITES  AND  BARROW 

G.  Weller  and  S.  Parker 
Routine  Monitoring  Program 

Radiation  fluxes  of  short-wave  incoming  (global),  short-wave  reflected 
and  net  all-wave  radiation  were  measured  continuously  through  sutler  1971  at 
both  Barrow  (27  May  -  5  September)  and  1-3  (22  April  -  29  October).  Epple, 
precision  pyranometers  and  Fritschen  net  radiometers  were  used.  Figure  1 
shows  the  installation  of  the  equipment  at  T-3  in  April  1971.  A  Davos  PD-1 
hemispherical  all-wave  sensor  can  be  seen  mounted  on  the  tripod  in  the 
foreground  of  the  photo.  This  instrument  was  used  on  occasions  to  obtain 
the  outgoing  all-wave  radiative  flux.  This  installation  was  replaced  by 
four  Eppley  pyranometers  (2  short-,  2  long-wave)  in  May  1972  and  has  been 
operating  since  at  T-3.  A  similar  set  of  sensors  was  used  on  the  AIDJEX  main 

camp  in  the  Beaufort  Sea  during  April  1972. 

Icing  problems  during  1971  were  severe  on  occasions,  more  so  at  T-3 
than  at  Barrow,  and  frequent  cleaning  of  the  Instruments  and  correction  of 
the  records  was  necessary.  During  1972,  this  was  made  unnecessary  by  the 
use  of  a  warm-air  blower  system,  especially  designed  and  built  for  the  pur¬ 
pose,  which  kept  the  sensors  completely  free  of  frost  deposits.  Table  1 
shows  a  typical  computer  printout  of  data  for  T-3  and  Barrow.  Monthly 
summaries  of  total  daily  integrated  radiation  values  will  also  be  listed 
in  the  final  tabulations.  Figure  2  shows  a  plot  of  typical  values  of  global 
and  net  radiation  when  skies  were  clear  at  both  stations.  This  occurred 
during  the  summer  solstice  on  21  June  1971  when  Barrow  (latitude  71"»  and 
T-3  (latitude  85”N)  were  approximately  2000  km  apart.  The  effect  of  latitude 
on  the  diurnal  variations  of  the  radiative  components  is  clearly  recognizable. 
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Installation  of  radiation  equipment  at  T-3  in  April  1971 
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These  basic  measurements  are  useful  in  synthesizing  a  radiation  climatology 
of  the  Arctic  Basin,  particularly  if  carried  out  over  longer  time  spans,  as 
proposed  by  us.  With  an  array  of  stations,  the  dynamic  nature  of  the  inter¬ 
action  between  cloud  advection  and  formation  and  the  radiation  balance  can 
also  be  studied.  Such  a  high-density  array  will  be  provided  by  the  full- 
scale  AIDJEX  station  grid  during  1974/75,  in  which  wa  propose  to  carry  out 
these  radiation  studies. 

Ice  Surface  Temperatures 

An  important  objective  of  the  field  studies  during  spring  1972  has  been 
to  extend  the  single-point  observations  at  the  stations  listed  above,  to  a 
greater  spatial  consideration  of  surface  inhomogeneities.  These  inhomogeneities 
are  present  in  the  form  of  different  types  of  ice,  snow  cover,  ice-free  and 
refrozen  leads  of  various  thickness,  melt-ponds,  hummocks  and  ice  ridges. 

Spatial  sampling  of  these  types  of  surfaces  have  included  measurements  of 
albedoes  and  radiative  surface  temperatures.  A  study  by  Maykut  and  Unter- 
steiner  (1969)  has  shown  how  even  small  differences  in  albedo  can  be  criti¬ 
cal  in  considerations  of  che  stability  of  a  sea  ice  cover. 

Radiative  surface  temperatures  were  measured  at  the  AIDJEX  camp  in  the 
spring  of  1972,  when  their  immediate  use  as  ground-truth  for  the  NASA  remote 
sensing  overflights  was  greatest.  The  basic  instrumentation  used  for  this 
purpose  was  a  PRT-5  infrared  thermometer.  Table  2  summarizes  the  results 
of  these  measurements.  It  is  clear  that  despite  the  relative  homogeneity 
of  multi-year  pack  ice  surfaces  in  spring,  large  temperature  differences 
can  exist.  This  means  that  single-point  observations  at  the  surface  have 
to  be  examined  carefully,  before  being  used  as  ground-truth  data  for  aerial 
remote-sensing. 
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TABLE  2 


Surface  Temperature  Differences  of  Multi-Year  Pack  Ice  in  Spring 


Weather  Conditions 

Surface  Temps  of  Sloping  Ice  Surfaces 
With  a  Flat  Surface  (°C) 

AT  max  (°C) 

Sunny  Side 

Shaded  Side 

Clear  sky,  calm 

+3.5 

-3.8 

7.3 

Clear  sky,  light  wind 

+1.2 

-1.8 

3.0 

Thin  clouds,  light  wind 

+0.7 

-0.3 

1.0 

0/cast,  strong  winds 

0 

0 

0 

This  tentative  table  represents  the  results  of  36  series  of  measurements; 
complete  analysis  of  additional  measurements  on  a  more  rigorous  basis  in  in 
progress. 

Radiative  Characteristics  of  Refreezing  Leads 

Much  more  dramatic,  of  course,  is  the  dynamic  temperature  change  of  re¬ 
freezing  leads.  This  was  measured  at  the  AIDJEX  main  camp  in  April  1972,  both 
on  naturally  occurring  and  artificially  created  leads,  and  is  shown  in  Figure 
3,  together  with  freezing  rates  and  albedoes  of  the  same  leads.  Mean  air 
temperatures  were  constant  at  around  -25°C  during  the  studies  so  that  steady- 
state  conditions  can  be  assumed.  Also,  n^  snow  fell,  or  was  blown  on  the 
refreezing  leads  during  the  period  of  observation.  Thermal  coupling  with  the 
water  under  the  ice  is  very  strong,  until  the  ice  reaches  a  thickness  of 
about  20  cm,  but  even  when  the  ice  is  a  meter  thick,  the  ice  surface  tempera¬ 
ture  is  still  up  to  2°C  warmer  than  the  near-surface  air  temperature.  A 
snow  cover  on  the  ice,  because  of  its  low  thermal  conductivity,  will  couple 
much  more  strongly  with  the  atmosphere,  so  that  temperature  differences 
between  air  and  lead  surface  will  be  substantially  reduced. 

Albedo  measurements  over  the  refreezing  leads  show  rapid  increases  with 
ice  thickness,  except  for  a  short  initial  decrease.  This  latter  effect  occurs 
when  the  mirror  reflectivity  of  the  open  water  surface  at  low  solar  elevations 
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RATE  OF  FREEZING,  SURFACE  TEMPERATURES  AND  ALBEDOS  OF  A  REFREEZING  LEAD 

AIDOEX,  APRIL  I97E 


is  reduced  due  to  ice  crystal  formation  on  that  surface.  From  these  values 
and  from  continuous  measurements  of  incoming  short-wave  and  long-wave  radia¬ 
tion,  radiation  balances  can  now  be  computed  for  the  lead.  These  are  shown 
by  four  days  following  freeze-up  in  Figure  4.  By  comparison  with  the  mea¬ 
sured  radiation  balance  of  multi-year  ice,  it  is  seen  that  the  open  lead 
looses  about  twenty  times  more  radiative  energy  to  the  atmosphere  than  Ln.e 
thick  multi-year  ice.  These  measured  values  agree  well  with  the  estimates 
by  both  Badgley  (1966)  and  Doronin  (1969)  for  early  April.  As  freezing  of 
the  lead  continues,  its  radiation  balance  approaches  that  of  the  multi-year 
ice  asymptotically.  Four  days  after  freeze-up  it  is  still  six  times  higher, 
the  multi-year  ice  balance  having  remained  constant  under  clear  skies  and 
identical  temperature  conditions.  The  importance  of  leads  in  the  heat  balance 
of  the  Arctic  atmosphere  is  thus  clearly  shown.  In  addition,  latent  heat 
fluxes  from  the  relatively  warm  lead  surfaces  also  contribute  heat  to  the 
atmosphere.  The  open  leads  may  well  supply  two  orders  of  magnitude  more 
thermal  energy  than  multi-year  ice,  as  suggested  by  Badgley  (1961). 

Light  Extinction  in  Snow.  Ice  and  Sea  Water 

Inhomogeneity  of  the  pack-ice  surface  also  strongly  affects  the  trans¬ 
mission  of  radiation  through  the  ice.  This  applies  in  particular  to  the  depth 
of  snow  on  the  ice  surface,  since  snow  has  a  much  higher  extinction  coefficient 
than  ice.  To  measure  the  extinction  processes  in  sea  ice  and  -ea  water, 
small  (2x1  cm)  selenium  photocells,  the  spectral  responses  of  which  cover 
the  solar  spectrum,  were  lowered  through  small  boreholes  in  the  typical  wind- 
packed  spring  snow  of  density  0.36  gm  cm  ^ ,  irregularly  covering  the  pack 
ice,  was  determined  by  placing  a  large  block  of  snow  over  a  photocell,  and 
reducing  its  thickness  by  cutting  sections  from  the  top.  Exponential  ex¬ 
tinction  occurred  in  all  three  media,  as  expected,  and  as  shown  elsewhere 
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FIG.  4: 


RADIATION  BALANCE  OF  A  REFREEZING  LEAD 
AIDJEX  CAMP,  APRIL  1972 
(CLEAR  SKIES,  TEMPS:  -  22®C  i  T>  -  30°C  ,  MEAN: -26°C) 


DAYS  AFTER  LEAD  FREEZING 


RADIATION  BALANCE  O  ©  ®  • 

ABSORBED  SOLmR  RADN.  □  D  [|  ■ 
NET  LONG -WAVE  LOSS  A  A  \  A 


MULTIYEAR  ICE 
(AIDJEX) 
BADGLEY,  LEAD 


by  Weller  and  Schwerdtfeger  (1967)  and  Weller  (1969). 


The  exponential  extinction  coefficients  for  snow,  Ice  and  water  were 
0.29,  0.012  and  0.0025  cm  \  respectively.  A  nomogram  can  be  constructed 
from  these  results  (Figure  5)  from  which  the  light  intensity  at  any  vertical 
level  below  the  surface  can  be  quickly  read,  for  any  combination  of  thick¬ 
nesses  of  snow,  ice  and  water.  For  example,  the  light  intensity  at  the 
ice-water  interface  under  10  cm  snow  and  2  m  ice,  is  0.7%  of  that  penetrating 
the  snow  surface;  for  no  snow  and  3  m  ice  it  is  about  3%. 

‘K 

Cloud  Albedoes 

The  arctic  radiation  regime  is  dominated  by  low  stratus  cloud  decks 
during  summer,  which  have  a  considerable  effect  on  the  radiation  balance 
of  the  surface  and  the  atmosphere  (Fletcher,  1965).  Our  efforts  to  define 
the  role  of  stratus  clouds  in  the  arctic  atmosphere  so  far  have  concentrated 
on  determining  the  physical  and  optical  properties  of  these  clouds,  through 
a  series  of  measurements  with  aircraft-mounted  sensors.  These  flights  took 
place  at  Barrow  with  light  aircraft  during  August  and  September  1971,  and 
again  in  April  and  May  1972,  flying  through  multi-layered  cloud  decks  up 
to  altitudes  of  4000  meters. 

Cloud  albedo  measurements  were  carried  out  at  the  same  time  as  cloud 
sampling  flights  (Section  III),  using  a  light-^weight  CSIR0  albedometer, 
mounted  forward  of  the  wing  and  attached  to  the  strut  of  a  Cessna  180  air¬ 
craft.  Figure  6  shows  typical  results,  indicating  the  wide  spectrum  of 
values  obtained  when  the  cloud  cover  is  thin.  Aircraft  flying  altitudes 
were  200  meters  above  the  cloud,  so  that  fifty  percent  of  the  sensor  re¬ 
sponse  comes  from  a  cloud  surface  area  400  meters  in  diameter.  This  is  a 
small  area  in  terms  of  resolution  elements  of  meteorological  satellites, 
so  that  spatial  variations  of  cloud  albedoes  due  to  holes  will  probably 
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S  5:  NOMOGRAM  FOR  THE  DETERMINATION  OF  LIGHT  INTENSITIES 
IN  A  THREE-LAYER  SYSTEM:  SNOW,  SEA  ICE,  SEA  WATER 

AIDJEX,  1972  ^8 

s  example  (heavy  line)  shows  light  intensities  for  10cm  of  snow  overlying  sea  ice  2  meters  thick. 


ALBEDO  (%) 

Figure  6.  Cloud  albedo  distribution  for  single  and  two-layered  cloud  deck 


not  be  seen,  except  in  the  mean  values.  On  the  other  hand,  satellites  such 
as  the  proposed  ERTS-series  will  have  several  resolution  elements  covering 
the  400  meter  diameter  area,  and  will  thus  see  the  holes.  The  measured 
values  appear  to  be  in  good  agreement  with  those  obtained  at  lower  latitudes 
and  reported  by  Fritz  (1957)  despite  the  fact  that  the  composition  of  the 
clouds  is  somewhat  different  (Section  III) .  It  is  clear  that  both  cloud 
thickness  and  the  presence  of  holes  in  the  cloud  will  affect  the  longwave 
'  as  well  as  the  shortwave  radiative  exchange  with  the  upper  atmosphere  and 
space.  The  latter  effect  is  discussed  in  Section  II. 

Ice  Surface  Albedos 

In  a  similar  fashion,  albedo  values  for  the  pack  ice  surface  were  deter¬ 
mined  at  the.  AIDJEX  camp  in  spring  1972.  The  instrumentation  in  this  case 
consisted  of  two  Eppley  precision  pyranometers  mounted  on  top  and  under  the 
rear  fuselage  of  a  Twin  Otter  aircraft.  The  flights,  primarily  photographic 
missions  at  regular  intervals,  consisted  of  one  run  of  10  km  length  at  12,000 
feet  over  the  AIDJEX  camp,  followed  by  j.0  runs  at  5,000  feet,  covering  an 
area  10  x  10  km  over  the  camp.  Finally,  a  low  run  at  1,500  feet  was  made 
over  the  camp.  As  with  the  cloud  albedo  measurements,  recording  was  on  a 
portable  single-channel  high-speed  line  recorder,  recording  incoming  radia¬ 
tion  for  5  minutes,  followed  by  10  minutes  of  reflected  radiation.  Mea¬ 
surements  were  only  made  when  the  sky  was  completely  clear.  All  data  are 
not  analyzed  yet,  but  initial  results  are  shown  in  Figure  7. 

As  far  as  these  tentative  results  allow  analysis,  the  following  state¬ 
ments  can  be  made.  The  mean  albedo  of  78%  for  multi-year  ice  agrees  well 
with  single-point  surface  measurements.  The  overall  frequency  distribution 
is  skewed,  due  to  the  presence  of  many  refrozen  leads,  generally  less  than 
a  meter  thick  and  more  than  1  km  wide.  From  the  histogram  for  all  observations 
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from  5,000  feet  height,  it  appears  that  they  cover  approximately  20%  of  the 
10  x  10  km  area  around  the  AIDJEX  camp.  Their  albedo  is  about  75%;  this 
relatively  high  value  is  due  to  snow  deposits  on  the  leads.  Smaller  fea¬ 
tures  with  different  albedos  are  not  seen  in  the  historgrams  -  also  the 
paucity  of  data  analyzed  so  far  does  not  yet  allow  comparison  of  data  from 
various  flight  altitudes. 

II.  EFFECTS  OF  ARCTIC  STRATVS  CLOUDS 
ON  THE  RADIATION  REGIME 

Gerd  Wendler 

The  arctic  summer  stratus  cloud  cover  is  one  of  the  most  persistent 
cloud  systems  of  the  world  and  consequently,  it  has  a  great  influence  on 
the  climate.  Its  importance  has  been  pointed  out  by  Fletcher  (1966)  and 
in  the  NAS  study  "Polar  Research"  (1970).  Theoretical  as  well  as  experi¬ 
mental  investigations  of  the  relation  between  incoming  solar  radiation  and 
clouds  have  been  carried  out  by  Albrecht  (1933),  Hewsen  (1943),  Neiburger 
(1948),  Haurwitz  (1948),  Vowinckel  and  Orvig  (1962a, b),  Loewe  (1963), 
Franceschini  (1968)  and  many  more. 

In  this  study  the  influence  of  stratus  clouds  on  the  incoming  radia¬ 
tion  was  investigated  lor  Barrow,  Alaska.  The  radiation  (incoming  and 
reflected)  was  measured  with  Eppley  pyranoraeters ,  while  cloud  observations 
were  obtained  from  the  Weather  Bureau  records  for  Barrow.  So  far,  the  study 
has  been  carried  out  only  for  the  month  of  June.  This  month  was  chosen 
because  (a)  the  solar  height  is  relatively  large,  so  that  the  energy 
received  from  the  sun  is  substantial  and  (b)  this  is  the  month  in  which 
a  great  change  in  the  surface  albedo  occurs.  The  experimental  site  was 
located  3  km  inland  from  Barrow  so  that  the  influence  of  the  Arctic  Ocean, 
which  can  have  a  very  different  albedo  than  the  tundra,  was  negligible. 

First  an  investigation  of  the  dependence  of  the  radiation  received 


on  ceiling  height  of  stratus  clouds  was  carried  out.  Figure  1  shows,  for 
three  classes  of  ceiling  heights,  no  systematic  difference  in  the  incoming 
radiation.  This  is  to  be  expected,  as  not  the  height  but  the  thickness 
of  the  cloud  cover  and  the  density  and  size  distribution  of  the  cloud 
droplets  will  determine  the  reflection,  absorption  and  transmission  of  the 
incoming  radiation. 

In  Figure  2.  the  amount  of  incoming  radiation  is  shown  as  a  function 
of  thickness.  (The  cloud  thickness  was  estimated  from  the  humidity  data 
of  the  radiosonde  ascents  at  Barrow.)  One  sees  that,  with  increasing  cloud 
thickness  and  especially  for  multi-layered  clouds,  the  amount  of  radiation 
received  on  the  ground  decreases;  however,  there  is  quite  a  substantial 


amount  of  scatter  in  the  data. 

Multiple  reflection  between  the  ground  surface  and  the  lower  portion 
of  the  cloud  is  of  great  interest,  since  it  has  considerable  influence  on 
the  amount  of  energy  received  at  the  surface.  If  one  assumes  that  a  cloud 
scatters  the  light  in  all  directions  equally,  which  is  approximately  true 
for  a  cloud  of  sufficient  thickness,  then  the  schematic  simplified  diagram 
given  in  Figure  3  is  applicable.  It  can  be  seen  that  the  energy,  which 
will  be  received  on  the  surface,  Isurface>  is: 


[  r  =  I  (1  -  a)  (1  -  a)  (1 
surface  u 


with  z  ■  (1  -  a)ct3 

2 

I  *  energy  received  above  the  cloud  per  cm  horizontal  surface 
a  ■  absorption  in  cloud 


u  ■  albedo  of  cloud 
3  ■  albedo  of  cloud 

One  can  see  that  the  multiple  reflection  becomes  increasingly  important 
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Fig  3  A  simplified  model  of  multiple  reflection  between  a  cloud 
layer  and  the  surface. 


with  increasing  cloud  and  surface  albedo  and  with  decreasing  absorption 
in  the  cloud. 

Calculations  for  the  incoming  radiation  were  carried  out  under  the 
following  simplifying  conditions: 

1)  The  incoming  radiation  on  a  cm  horizontal  surface  with  a  solar  elevati 
of  30°  was  estimated  to  be  0.80  cal  cm"2  min"1. 

2)  Tb'>  albedo  of  the  cloud  was  assumed  to  be  55%  (this  would  be  the  albedo 

as  seen  if  the  cloud  were  above  a  black  body). 

3)  The  absorption  in  the  cloud  was  assumed  to  be  7%. 

Under  these  simplifying  assumptions,  the  following  values  were  found 
for  three  days,  which  had,  according  to  the  radiosonde  ascents,  similar 
cloud  conditions: 


1.  Date 

2 

10 

29 

June  71 

2.  Hours 

11 

4 

5 

3.  Surface  albedo 
(measured) 

79 

58 

20 

% 

4.  Incoming  radiation 
(measured) 

.56 

.38 

.26 

i  -2 

cal  cm 

min  1 

5.  Mean  solar  angle 

31 

24 

21 

degree 

6.  *Incoming  radiation 
(measured) 

.54 

.47 

.37 

_2 

cal  cm 

min  1 

7.  *Incoming  radiation 
(calculated) 

*Lines  6  and  7  were  interpolated 

.55  .44 

to  30°  solar 

.36 

elevation. 

i  -2 

cal  cm 

min  1 

One  can  see  that  the  agreement  between  calculated  and  observed  values 
are  good.  The  decrease  in  incoming  radiation  of  32%,  which  was  observed 
for  similar  cloud  conditions  between  the  beginning  and  end  of  June,  can 
be  explained  solely  by  multiple  reflection.  At  the  beginning  of  June, 
the  ground  is  covered  with  dry  snow  and  the  albedo  is  therefore  high  (around 
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80%),  and  multiple  reflection  is  important.  On  June  10,  the  snow  cover  is 
melting,  which  lowered  the  reflectivity  by  about  20%,  and  decreases  the  im¬ 
portance  of  the  multiple  reflection.  At  the  end  of  June,  the  snow  cover  is 
melted  altogether  and  an  albedo  value  of  20%,  which  was  observed,  is  typical 

for  the  tundra  surface.  Now  the  multiple  reflection  is  unimportant,  which 

-2  -1 

explains  the  decrease  in  the  incoming  radiation  from  0.54  cal  cm  min 

-2  -1 

on  June  3  to  0.37  cal  cm  min  on  June  9.  Theoretically,  assuming  the 
cloud  conditions  were  identical,  the  decrease  in  the  incoming  radiation 
would  have  been  expected  to  be  even  slightly  larger  (35%) . 

Summarizing,  it  has  been  shown  by  an  example  that  multiple  reflection 
is  important  in  the  Arctic.  Therefore,  it  is  not  possible  to  give  typical 
values  of  incoming  radiation  for  a  specific  cloud  type  without  knowing 
the  conditions  of  the  underlying  surface.  If  one  would,  for  example,  cal¬ 
culate  the  radiation  balance  for  an  ice-free  Arctic  Ocean  by  using  the  ob¬ 
served  incoming  short  wave  radiation  data  obtained  from  ice  islands,  a 
serious  error  would  be  introduced. 

However,  more  work  has  to  be  done  in  this  direction.  So  far  observa¬ 
tional  data  of  only  4  days  were  used.  More  radiation  data  with  better  cloud 
observations  are  needed.  Also  the  model  is  too  simple  and  should  be  improved. 

III.  COMPOSITION  OF  ARCTIC  CLOUDS 
K.  Jayaweera  and  T  Ohtake 

Studies  of  the  composition  of  arctic  stratus  clouds  were  made  in  two 
series  of  experiments  at  Barrow,  Alaska;  the  first  series  was  performed 
in  September  1971  and  the  second  in  April  1972.  The  results  of  the  initial 
experiments  are  now  completed  while  the  analysis  is  in  progress  for  the 
second  series. 

The  cloud  particles  were  collected  continuously  in  flight  using  an 
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MRI  cloud  particle  sampler  extending  through  the  door  of  a  Cessna  180  air¬ 
craft  so  that  the  sampling  point  is  well  outside  the  disturbed  propeller 
wake  (Figure  1).  This  instrument  samples  air  at  the  rate  of  about  1.5 
i/sec  at  100-knot  flight  speeds.  A  4%  solution  of  formvar  in  chloroform 
continuously  coats  a  16  mm  clear  leader  film  just  before  it  exposes  itself 
to  the  cloud-air  through  a  small  rectangular  orifice.  Permanent  replicas 
of  the  cloud  particles  formed  in  the  hardening  formvar  are  later  analyzed 
in  the  laboratory.  We  found  film  speeds  of  4.5  inches  per  second  during 
sampling  the  most  satisfactory  rate. 

Clouds  during  September  were  rather  typical  of  the  summer  and  fall 
and  consist  entirely  of  stratus  decks  with  weak  or  no  updraughts.  They 
usually  are  multi-layered,  very  often  made  up  of  two  and  sometimes  three 
separate  layers.  Each  of  these  cloud  layers  are  about  1000  m  thick  and 
are  usually  separated  from  one  another  by  a  few  hundred  meters.  We  found 
clouds  extending  from  nearly  200  m  to  over  4000  ra  above  ground,  hence  we 
were  able  to  sample  clouds  at  temperature  levels  from  +2  to  -11C.  Sampling 
was  done  by  flying  at  a  fixed  height  or  constant  temperature  for  about  5 
minutes  so  that  about  500  liters  of  air  were  sampled. 

The  ice  crystal  concentrations  were  determined  by  analyzing  sections 
of  at  least  20  feet  of  film  at  a  time,  frame  by  frame,  using  a  stop  motion 
projector  and  counting  the  number  of  ice  crystals  within  each  section.  The 
ice  crystals  were  all  columnar  and  unrimed.  The  concentrations  of  ice 
crystals  at  -11  and  -7C  were  about  40  per  m3  and  10  per  m3,  respectively, 
and  at  -4C  it  was  less  than  2  per  m3.  These  values  are  comparable  to  the 
ice  nuclei  concentrations  at  the  corresponding  temperatures  measured  by 
collecting  air  samples  in  a  Millipore  filter  and  developing  them  using  the 
technique  by  Stevenson  (1968)  (Figure  2  and  Section  IV). 
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Concentration  per 


Figure  2. 


Comparison  of  Ice  Crystal  Concentrations  in  Arctic  Stratus 
and  Ice  Nuclei  Concentrations 
(At  Barrow,  Alaska  September  1971  and  April  I97£} 
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ope.  One  such  photograph  is  shown  in  Figure  3  The 
.  ,  s  e  3.  The  concentration  o £  the 

"P  aVeM8ed  ab°Ut  90  “3  «  levels,  it  must  he  mentioned 

-vsr.  that  the  concentration  v.rled  through  the  cloud,  owing  to  the  ' 

patchy  nature  of  the  latter 

er.  iron,  almost  zero  to  about  150  per  cm3.  The 

variation  ^t^dif fere  "  "  ’  ^  ^  dlStrlbUtl°"  “  ”«  =„„„  significant 

erent  l6VelS'  he"Ce  thc  °vera11  size  distribution  for  all 

d  ;V  8  18  Sh0Un  ln  -  course,  was  corrected  for  the 

fP  “  8Pr8ad  and  C°UeCtl0"  e£,1Cl^  <*  the  droplets.  The  3  „„  cutofc 
S  a  result  of  the  zero  collection  efficiency  of  the  sampler  for  droplets 
below  this  size.  P1  Cs 

IV.  ICE  NUCLEI  CONCENTRATIONS  IN  THE  ARCTIC  REGION 
T.  Ohtake  and  K.  Jayaweera 

Ice  nuclei  are  necessary  particles  to  form  ice  crystals  in  th 
Normally,  these  are  not  present  in  ab  H  atmosphere. 

nuclei  To  lnltl  t  1  ,UantltleS  3  *  condensation 

I.  TO  initiate  precipitation  in  cold  regions,  ice  crystals  should  be 

Present  in  supercooled  water  clouds,  even  in  small 

»  v-n  in  small  quantities,  l  ence  it  i8 

important  to  know  the  concentration  of  ice  nuclel  T 

.  .  nuclei.  In  arctic  regions  occa- 

onal  snowfalls  occur  even  in  sunder  and  further  research  into  global  h 
«or  the  influx  of  fresh  ice  nuclei  from  th  ‘  »cchanismt 

°  er  reelons  through  either  the  lower 

r  upper  atmosphere  is  needed  TKm  i 

neeoed.  The  importance  of  ic-  crv^fai  - 

crystal  aerosols  on  the 

arctic  radiation  regime  has  been  pointed  out  by  Fletcher  „966, . 

During  April  1972,  concentrations  of  ice  nuclei  were  observed  at  Barrow 
laska,  by  the  following  two  rn.tl.ods: 
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Figure  4.  Mass  and  number  distribution  of  Arctic 


1)  Cloud-settling  chamber  of  ice  nuclei  counting  (Ohtake,  1971).  Out¬ 
side  air  was  introduced  into  the  10  A  chamber  using  a  vacuum  pump.  After  the 
air  is  cooled  to  a  desired  temperature  (-15°C  or  -20°C) ,  supercooled  cloud 
droplets  are  introduced  into  the  chamber  to  give  water  saturation  exactly 
the  same  as  that  in  natural  clouds.  If  any  ice  nuclei  exist  in  the  chamber, 
these  will  changed  to  ice  crystals,  which  will  then  fall  into  a  supercooled 
sugar  solution  and  grow  to  visible  size.  The  number  counted  will  be  the  number 
of  ice  nuclei  in  the  10  ..iters  of  air .  Figure  1  shows  all  the  data  we  observed 
at  -15°C  and  -20°C;  they  are  compared  with  concentrations  of  ice  crystals  in 
the  atmosphere  (see  section  on  cloud  composition) . 

2)  Millipore  filter  method.  The  air  containing  ice  nuclei  was  drawn 
in  through  two  Millipore  filters  once  daily;  volumes  sampled  were  250  liters 
and  360  liters.  The  filters  ate  developed  in  a  manner  similar  to  that  of 
Stevenson  (1968)  in  a  processor  at  the  laboratory  in  the  Geophysical  Insti¬ 
tute,  as  follows: 

The  filters  are  cooled  to  -20°C  and  are  subjected  to  various  amounts 
of  moisture  in  the  processor  for  20  minutes.  The  number  of  ice  crystals 
formed  on  a  filter  indicates  the  number  of  nuclei  of  the  sample  volume.  This 
analysis  is  in  progress  now. 

A  postulate  for  a  transport  mechanism  for  ice  nuclei  will  possibly 
emerge  after  all  the  data  are  analyzed.  The  comparison  of  data  obtained  by 
the  cloud-settling  method  and  the  Millipore  filter  method  will  be  valuable, 
particularly  if  compared  further  with  the  measured  concentrations  of  ice 
crystals  and  ice  nuclei  in  the  undisturbed  polar  air,  because  of  complica¬ 
tions  in  the  measuring  techniques  concerned  with  ice  nuclei  concentiation 
(Bigg,  1970). 

Condensation  nuclei  concentrations  were  observed  in  addition  to  mea¬ 
suring  the  ice  nuclei  concentrations.  Condensation  nuclei  counts  consistently 
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Figure  1.  Concentrations  of  ice  nuclei  effective  at  -15°C  and  -20°C 
and  condensation  nuclei  (Aitken  nuclei)  t  Barrow,  Alaska. 
Note  that  the  concentrations  of  irp  nnrlo-i  at-  _onor> 


n/ce 

3000 


increased  when  winds  shifted  to  directions  from  inhabited  areas,  even  though 
absolute  numbers  in  unit  volume  of  air  were  quite  small  compared  with  other 
areas  in  Alaska.  On  the  other  hand,  ice  nuclei  concentrations  were  unaffected 
by  winds  from  inhabited  areas.  It  is  suggested  that  these  two  types  of  nuclei 
have  different  origins. 

Summarizing  the  results,  we  could  say  that  the  non-precipitating  stratus 
clouds  over  the  arctic,  within  thi^  temperature  range,  have  a  narrow  droplet 
spectrum  typical  of  strain-  clouds.  Ice  crystal  concentrations  are  low  and 
can  be  predicted  from  the  ice  nuclei  concentration. 

In  the  April  1972  experiments,  the  clouds  were  again  stratus,  except  on 
one  occasion,  when  strato-cumulus  was  observed.  On  two  occasions  the  sampling 
was  done  in  a  cloud  which  was  producing  snow  on  the  ground.  The  cloud  tempera¬ 
tures  were  lower  than  in  the  previous  experiments,  with  the  lowest  sampling 
temperature  being  -15C.  Preliminary  analysis  indicates  that  the  concentration 

of  water  drops  is  considerable  and  the  concentration  of  ice  crystals  is  higher 

3 

and  corresponds  to  that  expected  from  ice  nuclei  concentrations  (100  per  m 
at  -15C) .  In  these  temperature  regions  the  ice  crystals  are  plate-like  and 
dendritic.  Even  though  we  used  a  decelerator  to  reduce  the  air  speed  to 
one-half  the  aircraft  speed,  the  dendritic  crystals  were  all  found  to  have 
shattered  on  impact.  The  pieces  of  the  shattered  ice  crystals  were  found 
in  clusters  on  the  film.  However,  solid  hexagonal  plates  were  found  intact. 
Analyses  are  now  in  progress  to  determine  the  concentration  of  ice  crystals 
and  the  size  distribution  and  concentration  of  water  droplets. 
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V.  ABSORPTION  AND  SCATTERING  BY  ATMOSPHERIC  AEROSOLS 

G ,  Shaw 

Direct  solar  radiation  was  experimentally  monitored  at  several 
discreet  wavelength  intervals  in  the  visible  part  of  the  solar  spectrum, 
and  as  a  function  of  altitude  with  an  airborne  radiometer.  The  major 
scientific  objective  was  to  determine  the  height  and  wavelength  distribution 
of  the  optical  extinction  coefficient,  3(h,A),  arising  from  absorption  and 
scattering  from  atmospheric  aerosols  over  a  location  in  the  arctic  basin. 

Data  were  accumulated  near  Barrow,  Alauka,  to  an  altitude  of  5000  meters 
during  the  month  of  April  1972.  Supplementary  measurements  were  made  on  a 
Convair  990  aircraft  that  flew  north  from  Fairbanks  to  latitude  76°N.  These 
Convair  flights  extended  the  measuring  altitude  to  11  km.  This  report  summarizes 
initial  results;  the  theory,  instrumentation,  calibration,  and  experimental 
procedure  are  outlined  in  Appendix  I. 
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DISCUSSION  OF  DATA 


The  aerosol  extinction  coefficient,  3D(A,h),  was  obtained  by  the 
analysis  techniques  described  in  Appendix  I.  At  the  time  of  this  writing 
all  data  have  not  been  completely  analyzed;  however,  some  preliminary 
results  from  two  selected  flights  are  available  and  will  be  presented  in 
this  section. 

Flight  Number  4 

Figure  1  shows  the  analog  record  obtained  during  flight  number  4 
on  April  8,  1972,  extending  from  times  2025  to  2127  GMT.  The  sky  was 
completely  clear  ,  ground  temperature  was  -22.5°C  and  visibility  about 
20  miles.  The  flight  azimuth  was  held  at  110°,  climb  rate  was  adjusted 
to  500  ft/min  to  a  maximum  altitude  of  10,000  ft,  filter  at  5000  A  was 
used.  After  leveling  o££  at  10,000  ft,  the  aircraft  descended  at  a  rate 
of  400  ft/min  to  an  altitude,  of  approximately  50  ft;  during  the  descent 

O 

phase  the  7000  A  filter  was  used.  The  data  acquisition  ceased  at  a 
distance  of  85  miles  from  Barrow.  (Unfortunately,  data  could  not  be 
taken  during  the  return  flight  because  of  the  geometry  imposed  by  the 
aircraft  and  radiometer  part.)  The  solar  elevation  ranged  from  25.7° 
to  23.6°  during  this  flight. 

Figure  2  illustrates  the  aerosol  extinction  coefficient  at  A-5000 

O 

Angstroms  as  determined  by  the  analysis  method  described  in  Appendix  I 
section.  Plotted  adjacent  to  the  optical  extinction  coefficient  is  the 
temperature  profile.  The  radiative  inversion  extended  to  a  height  of 
2000  ft  with  an  average  lapse  rate  of  +34.0  °C/kms;  thereafter,  a  roughly 
isothermal  region  extends  to  about  6000  ft  elevation  after  which  the 
lapse  rate  approached  the  moist  adiabatic  value  of  -6°C/km  for  higher 
altitudes.  There  is  some  slight  structure  in  the  temperature  curve 
with  a  relative  maximum  (warm  region)  near  3000  ft.  The  optical 
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FLIGHT  NO.  4 
BARROW,  ALASKA 
APRIL  8  >  1972 
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Figure  1.  .  Analog  record  of  radiometer  output  voltage  (proportional  to  the 
aolar  intensity  at  5000  A)  obtained  during  flight  number  4. 
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extinction  coefficient,  3(h,5000  A),  shows  a  general  decrease  with  alti¬ 
tude.  The  individual  points  have  rather  considerable  scatter  which  are 
mainly  carried  by  error  in  reading  the  slope  of  the  I(X,h)  vs  h  curve. 

The  tendency  for  the  extinction  coefficient  to  peak  near  1500  ft  altitude 
seems  to  be  real;  other  short  periodicities  are  probably  noise  and  can  fa- 
averaged  out.  A  line  drawn  through  the  data  points  gives  an  approximate 

O 

empirical  expression  for  3(h,5000  A)  as 

li 

3(h)  =  0.410  e  H  (1) 

H  =  aerosol  scale  height  =  1.50  km 

The  total  optical  depth  (at  5000  A)  can  be  obtained  from  eq.  4  (Appendix 
I)  and,  at  the  maximum  altitude  of  10,000  ft  is  found  to  have  a  numerical 
value  of  0.231.  Subtraction  of  «-he  optical  depth  arising  from  mole¬ 
cular  scattering  above  10,000  ft  (tr=0.099)  and  optical  depth  arising 
from  ozone  absorption  (tg=0.014)  yields  an  optical  depth  of  0.117  above 

10,000  ft  elevation. 

Atmospheric  aerosols  are  distributed  in  two  height  regions,  namely 
in  the  lower  troposphere  and  also  within  the  stratosphere.  The  aero¬ 
sols  in  the  stratosphere  may  give  rise  to  the  majority  of  the  extinction 
above  10,000  ft.  To  estimate  the  magnitude  of  the  extinction  caused 
by  the  stratospheric  aerosols  we  assume  that  the  troposphere  extinction 
coefficient  can  be  estimated  for  heights  greater  than  10  kft  by  extra¬ 
polating  the  exponential  relationship  given  by  eq.  1.  Thus,  the  total 
aerosol  depth  above  10,000  ft  (2.94  km)  arising  from  the  exponential 
decrease  of  tropospheric  aerosol  is 

t (h>10  kft)  =  0.41  e  h/H  oh  =  0.092  (2) 

z.94 

Therefore, , the  residual  optical  depth  attributable  to  extinction  by 
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the  stratosphere  aerosol  would  be  t  »  0.117  =  0.092  «*  0  025 

strat 

Admittedly,  the  above  method  for  estimating  the  extinction  arising 
from  stratospheric  aerosol  is,  at  best,  somewhat  crude  and  almost  cer¬ 
tainly  the  height  profile  of  the  optical  extinction  coefficient  is  com¬ 
plex  and  cannot  be  well  represented  by  extrapolation  of  the  lower  level 
extinction  coefficient.  However,  the  method  is  plausible  and  offers  a 
rough  way  to  deduce  the  stratospheric  extinction. 

The  most  noteworthy  point  about  the  atmospheric  optical  parameter 
as  inferred  by  the  data  accumulated  on  flight  4  is  the  unexpected 
high  values  for  atmospheric  turbidity.  The  data  yields  a  total  value 
of  optical  depth  equal  to  .670  for  a  wavelength  of  5000  A  and  at  ground 
level.  After  subtracting  out  Rayleigh  and  ozone  absorption  one  is  left 
with  td=.510.  This  value  is  to  be  compared  with  a  standard  mid-latitude 
atmosphere  tabulated  by  Elterman  Q968)  who  gives  a  value  of  1^=0.264. 
Moreover,  the  high  value  of  optical  depth  found  at  Barrow  is  completely 
at  odds  with  actinometric  data  acquired  during  the  spring  and  summer 
of  1970  and  1971  at  McCall  Glacier  (69°18’N,  143°48'W)  in  northeast 

O 

Alaska  where  an  average  value  of  Angstrom  turbidity  coefficient  was 
found  to  be  approximately  0.020  (Shaw  and  Wendler,  1972)  which  corre¬ 
sponds  to  an  optical  depth  due  to  dust  at  A=5000  A  of  0.049  (assuming 
-1.3 

a  A  wavelength  dependency) .  At  an  altitude  of  1700  meters  corre¬ 
sponding  to  McCall  Glacier  our  measured  aerosol  optical  depth,  t^, 
was  0.220  or  a  factor  of  4.5  larger  than  the  background  level  as  re¬ 
corded  at  McCall. 

Inspection  of  other  data  acquired  during  April  shows  that  the  op¬ 
tical  depth  obtained  during  flight  4  was  unusually  large.  The  reason 
for  this  particularly  large  value  of  aerosol  optical  depth  is  unknown. 
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However,  some  clue  to  the  cause  may  be  found  by  Inspecting  the  synoptic 
m?ps.  Figure  3  shows  the  synoptic  situation  at  the  surface,  at  the 
850  mb  level  and  at  the  500  mb  level  during  flight  number  4.  By  in¬ 
specting  previous  maps  and  by  comparing  the  temperature  profile  trend, 
it  was  determined  that  a  cold  polar  air  mass  had  recently  been  in¬ 
truded  by  a  relatively  warm  continental  southerly  flow.  In  particular, 
the  circulation  near  the  surface  v^s  from  160°  azimuth  from  interior 
Alaska.  At  higher  levels  the  winds  veered  to  clockwise  to  290°  at 
850  mb,  and  about  10°  at  500  mb.  It  is  conceivable  that  the  high  tur¬ 
bidity  had  its  origin,  at  least  in  part,  because  of  the  iow  southerly 
flow  paths  over  interior  Alaska.  During  the  period  of  measurement  the 
ground  was  exposed  to  the  south  of  the  Brooks  Range  in  northern  Alaska 
and,  hence,  it  is  possible  that  the  large  dust  content  could  have 
originated  from  interior  and  southern  Alaska. 

The  high  turbidity  during  flight  4  could  also  be  partly  due  to 
sub-visible  cirrus  ice  clouds  associated  with  the  strong  low  pressure 
system  centered  at  some  400  miles  to  the  southwest  of  Barrow.  As  the 
flight  progressed  toward  a  bearing  of  80°  true  the  airplane  would 
move  further  away  from  the  low  pressure  system  and  thus  may  encounter 
clearer  air.  This  seemed  to  be  the  case  as  indicated  by  a  record 
of  6olar  intensity  obtained  while  flying  at  constant  elevation.  Thus, 
the  large  apparent  turbidity  may  be  due  partly  to  ice  crystals  in  the 
thin  cirrus  clouds. 

Flight  Number  5 

Flight  number  5  occurred  on  April  10,  1972  from  1422  to  1553  AST. 
Sky  conditions  were  clear;  visibility  approximately  50  miles;  wind  ENE 
at  15  knots;  ground  temperature  -17°C.  After  gaining  altitude  and 
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Figure  3.  Synoptic  situation  at  (a)  surface;  (b)  850  rob  level;  and  (c)  500  rob  level,  for  the  vicinity 
of  Barrow  during  the  time  period  of  flight  number  4. 


carefully  scanning  the  sky  some  very  tenuous  cirrus  were  visible  near 
the  azimuth  of  the  sun  but  considerably  below  the  sun.  The  solar  ele¬ 
vation  angle  during  the  flight  varied  from  23.5°  to  19°. 

O  O 

Optical  data  were  accumulated  at  5000  A  and  at  7000  A.  Evidence 
of  horizontal  inhomogeneities  encountered  during  the  acquisition  of 

O 

the  data  at  7000  A  render  its  value  somewhat  questionable.  However, 

O 

the  data  acquired  at  5000  A  during  both  ascent  and  descent  seem  to  be 
consistent  and  thus  these  data  will  be  discussed. 

O 

The  extinction  coefficient  arising  from  aerosols  at  5000  A  is 
shown  in  Figure  4.  Although  there  is  a  considerable 
scatter  in  the  individual  data  points,  the  trend  of 

decreasing  volume  extinction  with  increasing  altitude  is  quite  apparent. 
A  straight  line  fit  through  the  data  yields  a  scale  height  of  1.1  km 
which  is  somewhat  smaller  than  the  scale  height  found  one  day  earlier 
during  flight  number  4. 

O 

The  aerosol  optical  depth  evaluated  at  5000  Angstroms  at  the 
ground  and  at  peak  altitude  are  as  follows 
ip  ■  0.280  at  ground  level 
Tp  ■  0.062  above  12,000  ft  (3,52  km) 

The  values  above  are  about  half  of  the  corresponding  values  found  the 
previous  day  during  flight  number  4.  There  is  some  tendency  for  the 
aerosols  to  be  trapped  where  the  temperature  structure  is  concave  to 
the  right;  however,  the  inherent  scatter  of  the  individual  points  makes 
any  definite  inferrence  somewhat  difficult.  Again,  the  scatter  of  the 
points  is  mostly  caused  from  errors  in  reading  dl/dh  from  the  analog 
record  -  the  process  is  statistically  unbiased,  however,  and  the  average 
trends  are  meaningful.  The  large  variance  in  individual  points  is 
basically  due  to  noise  added  by  the  process  of  differentiation. 
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AEROSOL  EXTINCTION  COEFFIENT  (KM""1) 


Optical  extinction  coefficient,  8(5000, h)  obtained  during  flight 
number  5  on  April  10,  1972  near  Barrow,  Alaska  R 


The  residual  optical  depth  above  12,000  ft  0*52  km)  is  probably 

at  least  partially  caused  by  the  optical  extinction  by  the  stratospheric 

aerosols.  Assuming  the  continuation  of  the  exponential  dependency  of 

8(h)  to  heights  above  3.5  km  would  yield  an  optical  depth  equal  to 

0.012  above  3.5  km.  Thus,  it  would  be  inferred  that  td  stratosphere  - 

0.062  -  0.012  **  0.050.  It  becomes  very  difficult  to  accurately  assess 

optical  depths  of  this  small  magnitude  and,  in  fact,  an  error  analysis 

shows  that  the  residual  stratospheric  extinction  is  subject  to  an  error 

of  the  order  of  t  ■  0.050  +  0.025. 

strat 

Convair  990  Flights 

Data  from  the  Convair  990  flights  have  not  at  the  time  of  this 
writing  been  subjected  to  analysis.  Briefly,  two  flights  were  made, 
both  originating  from  Eielson  Air  Force  Base  and  flying  to  a  grid 
near  the  AIDJEX  campsite.  The  flights  extended  to  76°N  latitude  and 
encompassed  an  altitude  range  from  40  ft  ms£  to  37,000  ft  (10.8  km). 

O  P 

The  optical  measurements  were  made  at  7000  A  and  at  5000  A  through 
one  inch  thick  plexiglass  ports  with  vertical  angle  of  +14°.  Unfor¬ 
tunately,  a  great  deal  of  flying  time  was  spent  either  flying  almost 
into  or  out  of  the  solar  azimuth  and,  since  data  could  only  be  accumulated 
at  times  when  the  sun  was  at  nearly  right  to  the  aircraft  bearing,  only 
a  limited  amount  of  data  were  accumulated.  However,  these  measurements 
that  were  obtained  upon  cursory  inspection  seem  to  be  of  good  quality 
and  most  likely  will  provide  more  useful  information,  especially  in  re¬ 
gards  to  the  stratospheric  aerosol  extinction. 

SUMMARY 

The  airborne  method  of  studying  the  vertical  distribution  of  at¬ 
mospheric  aerosols  has  provided  a  considerable  amount  of  valuable  in¬ 
formation.  The  optical  extinction  coefficient  is  found  to  vary  in  the 
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S.  A,  Bowling 

Introduction 

Some  work  was  carried  out  during  the  year  „„  several  th  '  . 

Of  rediative  transfer  in  the  arctic  at  ,  a8PeCte 

previously  developed  for  ice  fog  ^  Infrared  Cra"Sfer 

a  clear  air  la  J  "aS  exPa"<‘ed  to  allow  for 

ar  air  layer  below  a  clo-J  deck  in  +  * 

in  the  r  dl  lon*  improvements  were  made 

in  the  temperature-dependence  of  the  model  UsW  a  , 
size  di  *■  8  8  m0del  and  a  measured 

distribution  for  arctic  stratus  clouds  (Section  III,  «  h 

ta*  ~  Without  the  Cloud  deck  were  c  „ 

were  computed.  In  addition,  an 
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approach  has  been  developed  which  allows  computation  of  fluxes  in  an  atmo¬ 
sphere  with  several  cloud  layers.  Programming  on  this  model  is  in  progress. 
Finally,  some  rough  calculations  have  been  made  on  the  effect  of  a  thin,  purely 
scattering  cloud  layer  on  the  total  energy  absorption  of  the  earth-atmosphere 
system  at  various  latitudes  and  seasons. 

3-Layer  Model 

The  3-layer  model  as  it  currently  exists  allows  for  two  clear  air 
layers  with  an  absorbing  and  scattering  layer  sandwiched  between  them.  Car¬ 
bon  dioxide  and  water  vapor  are  accounted  for  throughout  the  29-point  sounding, 
and  ozone  is  considered  to  provide  an  upper-boundary  input  in  the  10  ym 
region.  A  temperature  discontinuity  at  the  ground  is  allowed,  as  is  a 
wavelength-dependent  ground  albedo,  but  neither  of  these  features  has  been 
utilized  to  date.  Cloud  absorption  and  scattering  cross-sections,  as  well 
as  generalized  absorption  coefficients  for  water  vapor  and  carbon  dioxide, 
are  specified  over  63  wavenumber  intervals.  Interval  widths  are  40  cm  1  over 
most  of  the  spectrum,  with  10  cm  1  resolution  over  the  CO2  band  and  200 
cm"1  bandwidth  on  the  shortwave  side  of  1600  cm"1.  Variation  of  gaseous 
absorption  parameters  with  temperature  is  handled  by  Elsasser’s  (1960) 
parabolic  correction,  assuming  the  temperature  is  the  mean  temperature 
of  the  total  absorption  path  under  consideration.  The  absolute  accuracy 
is  probably  comparable  to  that  of  Elsasser's  tables  as  corrected  by 
Zdunkowski  et  al.  (1966);  the  difference  between  calculations  with  and 
without  clouds  is  probably  better  than  this.  Calculation  of  fluxes  and 
cooling  rates  over  the  29  sounding  intervals  requires  less  than  5  minutes 
of  CDC  7600  computer  time. 

A  standard  Mie  scattering  coefficient  program  available  at  the  National 
Center  for  Atmosphe  'r  Research  has  been  modified  to  accept  complex  refrac¬ 
tive  indices  for  the  63  wave  numbers  used  in  the  radiation  program  and  any 
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measured  or  theoretical  size  distribution,  to  calculate  the  absorption  an<< 
flux  backscatter  coefficients,  and  to  punch  the  results  on  cards  in  a  format 
suitable  for  direct  inclusion  in  the  radiative  transfer  program.  A  measured 
size  distribution  for  arctic  stratus  clouds  (Section  III)  was  combined  with 
the  complex  refractive  index  for  water  given  by  Irvine  and  Pollack  (1968) 
in  this  Mie  program.  The  results  are  shown  in  Figure  1.  An  August  sounding 
from  Point  Barrow  was  then  used  to  calculate  cooling  rates  with  and  without 
this  cloud,  the  cloud  height  and  depth  being  determined  by  the  saturated 
layers  in  the  original  sounding.  Latent  heat  effects  were  included  in  cal¬ 
culating  the  coo  Ling  rates.  Results  are  shown  in  Tables  1-3.  Flux  values 
are  in  cal/ (cm2  12  hours)  and  cooling  rates  in  °C/12  hours. 

In  general,  the  results  confirm  a  tendency  noted  in  every  computation 
made  thus  far.  The  rate  of  energy  loss  decreases  at  the  ground,  from  the 
atmosphere  below  the  cloud,  and  from  the  earth-atmcsphere  system  as  a  whole 
(assuming  the  cloud  temperature  is  below  that  of  the  ground).  Cooling  rates 
above  the  cloud  are  enhanced,  and  the  cloud  itself  may  show  increased  cooling 
rates  (for  low  clouds  and/or  low  surface  temperatures)  or  heating  (for  high, 
cold  clouds  over  warm  surfaces).  The  first  case  is  probably  always  true 
for  arctic  stratus. 

Visible  Radiation 

A  very  simple  model  of  the  effects  on  solar  radiation  of  a  thin  cloud 
was  developed.  The  model  included  multiple  reflections  between  cloud  and 
surface,  and  integrated  polar  angle  over  a  24-hour  day.  Cloud  particles 
were  assumed  to  be  pure  scatterers  and  a  single  scattering  event  (from  a 
cloud  particle  or  the  surface  of  the  earth)  was  assumed  to  destroy  all 
directional  characteristics  of  the  incoming  radiation.  The  albedo  of  the 
underlying  surface  was  assumed  independent  of  solar  zenith  angle.  The  model 


WAVE  LENGTH  [  p. m) 


dflrt  CflA  Qflft  innn  iAj>. 


FLUXES  AND  HEATING  RATES 


J 

I 


I 


W  N.  K  FO'Hlft(\j0>'U«£,K.thC)cO  4=3 
aeoNHQOJ(M(\j»oM*5tOk£!Hi 
pjj-CTflOJ-u'N.®trirvcri4''Hooir> 

•  *  •  i  •  ♦  •  k  •  i  h  i  ]  r 


i  I  !  !  I  •  '  / 

iOITicO  Jt  THs}'(TiO(MrOer5LlSfO»-:cr».* 
Q  ^  C\J  (\J  HHHNrtWHtVIQ'^OO 

OOOOQOUOOOoaOQOO 


rtci 

o  o  □  o 

•  •  •  ]• 

OOBIH 

O  O  O'  LA 
■H  C\J  l\J  LA 


O  O  O 
h  t»-  K 


•  »  •  > 

ST  IT  £  x: 
aoao 

•  •  •  < 

oootp 
a  o  O' 
*}«  Al  Al 


r  r  £  r  r  r 
i  a  a  O  O  O  O 
•  ••*•• 

I  LA  A»  cO 

r»  eo  o  (vj  vo 

<Olft  WIPOH 
r-L  CU  J-  .}■  LA  LA 


Al  f 

i 

l 

i 

O  O  o  I 


•  »  • 

r  £  £ 

o  □  a 
•  •  • 

LA  IA  O' 
IAN.N. 
CO  LA  A. 
>0  A.  rv 


>000600 

1  k"  >r*  •—  •“  •—  ►» 

.  !.  .  L  .  1 
£  £  £  £  £  £ 
OOOOOO 
•  ••••• 
H  IA  J’HJ’ 

IAMOON4 

O  CO  1A  K>  O'  O 
•H  r*  Al  Jh  LA 


OOO 


•  •  • 

£  £  £ 
OOO 
•  •  • 
CO  LA  IA 
vOIAN 
r*  CO  LA 
LA  i0  A. 


A.Mo^J-OlArOlA-Hj-iAO'K.otij 

OlflO'O'  ■tfODNWO'J'QJ-J-OM 

J-cOLAlALA^-O'cO'H'HfOLAlAO'CA-H 

.^iJI!***,**l,**»**i' 

lAvD(T>lA  OHv0O(\jO(\JtAO(J>(^O 
©«(00>ooHr0iA0  0^k^)a3  4- 
;  I  HrLHHHHHHHHHW 


N^HtOibrOhliDlAo^ocPlAoM# 

o^Aiwroty  -h  1  ^swio^sorirt 

iTl  1  !  |  hnnnw^  j 


i  _  .1  .  1  1  1  1  1  1  I 

ooooQooaooooooaoo 

•  OOO^HlOHIAAitH>t(OLAlAO>0 

£  j  000'1/\®U\Aioon^iDIAAAo 

*  I  ■HAJAJLAO'OeoiAfOO'O-HcOiAA.LA 

N  !  I  T*  -H  AJ  J-  *  LA  LA  V0  N  A.  l0 


29 


August  11,  1970,  1200Z,  No  Clouds 


FLUXES  ANO  HEATING  RATES 


I 

I 

I 

I 

I 

I 

!o 


I  “ 

K 

ix 

i 

i 

i 

i 

!> 

is 

x 
3 

U. 


to 

h- 

M 

r 


or 

w 


-T 
O' 
X  co 
3  • 
_l  J- 
U.  I 


N. 
O  -H 


l 

i 

l 

l 

i 

i  • 

|X 

i  • 

jN 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 

i 


1 

1 

1 

1 

1 

1 

!tn 

!  1 

1  1 

1  1 

1  1 

1  1 

1  1 

4  <\J  ^'ITN  vJDVl  C 

1  1 

!  1  1 

1  1 

1  1 

1  ! 

1  1 

oj  <n'M  in; 4-  * 

1 

1 

1 

co  o' 

N.  OJ  40  O  in  o  fO  J 

o>  o'  s  ro  co  \C 

'•H  J-! 

! -h  o,  ro  in  o'. co  <c 

in  id  O'  >T!«-i  a 

>••!••!•« 

0  in; 

• 

•  •! 

:  i 

•  ii  ^  *  » 

l  1  j  1  1 

i 

i 

i 

i 

i 

J 

j  t  j 

i  i 

i  1 

1  1 

1  1 

1 1 

i 

i 

i 

lift 

i  i 

i  i 

i  * 

IHvOlM  CT»1  <M  O 

1  1 

1  1 

1  1 

cvj  ro'k  in'ro  h 

O  4*1 

ia  c 

■h  cvj  <nj  nj  nj  f'j 

H  H  H  PJ  O  O 

O  O' 

Jo  c 

a  a  o  nj  o  o 

O  CJ  Ci  0  0  0 

0  oj 

1  • 

1  V  #|  •  V,  •  • 

•  •  ! 

!  I 

i 

i 

i 

i 

i 

J  Hj 

1  1 

1  1 

1  1 

i  i 

:  • 

1  1 

1  1 

I  1  1 

1  ! 

i 

i 

I  •  • 

Jr  r 

i  i 

t  i 

i  i 

r  r  r  r  r  r 

i  1  1 

1  1 

1  1 

t  t 1  *  t1  «  t 

r  r  'r  x'x  x 

•  •! 
r  rj 

io  a 

o  o  o  o  o  o 

000000 

O  Oi 

••••!•• 

•  •  ■ 

Jo  o 

O'  H  CO  *4  in  N. 

j’  h  j-  co  in  in 

O'  0! 

io  o 

o'  in  e  mN  ®  a  w  Jt  veins-  n.  oi 

!*h  <\j 

iW  in,  o'  acoinroo'o*-ia3  inn.  in! 

1 

i  HiH  OvJ 

■t  ^  in  in  id  s 

s  io> 

1 

1 

i  l 

i  I 

1  1 

1  1 

1  1 

mJ 

1 

!o  o 

i  i 

o  oo  OO  o 

1  1 

0  o'o  o!o  0 

1 

0  oi 

jt-  h- 

1 

t-jH-  h- 

1  1 

t-  t-  H-  t-  t-  t- 

h-  h~i 

1 

!  •  • 

tr  r 

x  x  x  rt  r 

•  •  1  •  •  1  •  • 

x  x.  :x  r  x  x 

•  •! 

r  r! 

]o  a 

o  o  o  oo op o  o  o  o  o 

0  o] 

.  . 

a  a 

o  O'  H  O  rlin 

n.  j’  *4  j-  to  in  in  o' 1 

a 

o  cr-  in  co  in  too  nj  4  io  in  n  s 

<H 

nr  nj  in  a'  o  <0 

in  ro  o'  0  *4  to 

in  si 

i  M  ^ 

1  1 

1  1 

<u  j<*  in  ;tn  to 

1  1 

1  1 

N.  N.J 

1 

l 

1 

O'  vO 

1  1 

1  1 

1  « 

(VJO'IrMAWyD 

1  1 

1  1 

1  1 

(TkinH^o^ 

i 

1 

1 

in  co ! 

<M  vO 

«x  «-i  in  co  ro  |v0  *-r\0  (\j  ro  <vj  n.  1 

J-  ro 

NJO^iN  J- 

to  O'  O  CJ  (\J  iD;in  S1 

•  • 

•  •••:••! 

•  •  •  •  •  •  1 

#  •  1 

in  in 

•r  01  n  o'  s  ^ 

NOIOJJHO1 

0  oi 

1  • 

1  M  «;on 

-tinininvDvD'ioro 

Jt-I  H 

1  1 

1  1 

1  1 

^  H  iH  H  iH  H 

1  1 

l  l 

1  l 

H  OJ] 

1 

1 

J-  <H 

1  1 

1  1 

^  1  1 

40  vO  fO  k  ivD  in 

l  l 

1  1 

1  l 

0  rl«0  (T  in  O 

1 

l 

1 

*>  <n  | 

•  • 

•  •!  •  ••  •  • 

•  •  •  •  J  •  • 

•  •  1 

m  ro 

re  nil  -Hit 

s  ro  to  -y  s-  0  '< 

h  ro] 

1  !  1 

h  pj  ,nj  nj  ro 

j-  J’i 

1  1 

1  1 

1  1 

1  1 

1  1 

1  I  1  1  III  | 

1  1 

1  1 

1  1 

1  1 ; 

1 

1 

1 

o  o 

1  1 

1  1 

1  1 

0  0  Jo  0  Jo  O 

1  1 

1  1  1 

1  1 

3  0  p  0  0  0  < 

1 

1 

1 

1 

3  O' 

•  • 

•  •  «i 

•  • ! 

o  o 

o'  <h  co  «-t  in  n> !. 

?  *4  90  in  in  <y»  o. 

o  o 

cn  in  ®  in  s  ©  orjj-vcms  r-  a; 

H  0  J 

m  in  o'  a  «o  in  ro  o'  0  *-t  c©  in  s.  in  1 

1  H  -H  0J  • 

t  ^  in  in  no  s  r 

-  vo; 

1  ^  1 

1  1 

(M  1 

+75  '*’! 

1  1 

1 

1 

"J 

1  1 

1 

d  | 

1  1 

1 

-J  ! 

(  1 

1 

U  ! 

1  1 

1  1 

1  l 

1  1 

1  1 

1  1 

1  1 
;  1 

54 

1 

1 

1 

1 

1 

1 

co 

B 


o 

10 

a 

o 

u 

T3 

O 


4J 

(0 

3> 

is 


M-UXES  AND  HEATING  RATES 


i  UN  O  {\j  <M 

|  K.  <\l  ®  O 
HO.fOlfl; 


;ooooofo;oooo;oo|SSSo! 

1 1  !  i 

!  i  i 

I*.*;..:.  I  i 

n  n  S  r!r  r  r  r'r  E'E>  t 

jOO  OCT  a  0,00,0  ojo  e.oo.o  oj 
!°  ®  “Jl®  UVK  ®  o  pj ij-  un  r.  o 

I*  n  ,vj  *  ®.®  £  w  T  °  *  *'*-  “i 

1  -t  4'iUN  IP  h.  iO I 


P  ™  ■  s  W 

jtltl 


i 

il  i 

I  i  j  • 

un;  «  oVpj  ov^  uv  4-  «o  oi 

S  !2  f?  W;«0  u>  H  J-! 

^  fO,  ®  SO  Lf\  t£)'  0>  H  O  O  UN,1 

•*  **1  .*  .*  •  •*••!•••  •! 

•  •  vfl,  *! 

•  ii  i 


i  i  1  i 

un  cm  *-i  o'rj  ro'K  uVfo 
o  un  pj  PJ,^  h'h  rvlS 


jO  O  O  o! o  oo  o  o  o'o  o 

!*“*“•“  *”1'”  •— *->-  *-! 


O  a  o  o 


|3C  X  JC  X[2Z  X|X  y  |  v  y  Jr  r'r  y'^- 
|o  o  o  O  o  O  O  o  5  O  S  S  a  o  S  £ 

j®  gig  £  £  «  ^  uvp!  ^‘4  ..*«  p 

i  ^)PJ  <\j;un  m,S  e”^Si  oiS  Sms 


JH  "HjPJ  4'  !«4  UN 


UN  iO!N.  K. 


LSiss^^^s^sisi’gs'sfc 

"i”  2°.  3  ~  ^  ® 

Cl  f?  T  "f  CJ  't  <Jl  N-  ^'.W  H  ro  loiri  o  o  o 

“*  '  !  '  1  '  1  «0|O  PJjLf  UN  UN  UN;*C  UJ  vO  Si 

!  j  ^  HM  Hiw  r);r(  rIM  (Mj 

!  !  !  '  !  ! 


NjLf  *4  «  vojro  k!«  un'o  h'«  ONjUN  O  fo  o>! 


JO  ^JPJ  ro  ro  pj,  ^ 

I  *  I  I 


I  t  ro  ro  *T:N.  oj«4  FC  I 
!  »  k*  y  y  ytr?  4-1 


0:0  O  O  o' 


|  w,  WWW  •  •  I  #  ~  ~| 

H  4-  0'iAiA  m'a! 
j*  |2  5  ®  °  W  ®  UN  N.  JN.  oj 

V?  '  nT2:2S;i2T°^i®un'l unI 

I  J  ’H;’H  pjj4-  -4  UN  UN  vD  P- 1 N.  U>! 


31 


TABLE  3 

August  11,  1970,  1200Z,  Cloud  Density:  80  drops/cm' 


was  run  for  ground  albedos  from  .05  to  .85,  10°  latitude  intervals,  and 

cloud  optical  depths  from  .01  (thin  haze)  to  1  (thin  cirrostratus)  at  both 

solstices  and  the  equinox.  The  results  indicated  thr.t  for  an  optical  depth 

of  1,  the  absorbed  shortwave  radiation  would  be  reduced  by  almost  300  cal/cm 

2 

day  at  20°  in  the  summer  hemisphere  and  by  150  cal/cm  day  at  the  summer  pole. 

The  latitudinal  albedos  of  Vonder  Haar  et  al.  (1971)  were  used.  If  the  summer 

pole  albedo  (witho  the  thin  cloud)  were  reduced  from  .60  to  the  global 

average  value  of  .30,  the  reduction  in  absorbed  radiation  could  be  as  high 
2 

as  350  cal/cm  day.  The  infrared  effect  of  such  a  cloud  will  depend  criti¬ 
cally  on  cloud  height,  but  calculations  with  cirrus  at  the  tropopauae  and 
at  the  -40°  C  isotherm  suggest  much  less  contrast  between  summer  and  winter 
hemispheres  in  the  longwave  component.  The  overall  effect  of  a  thin  cloud 
layer  would  appear  to  be  a  reduction  in  the  pole  to  pole  energy  budget 
gradient,  at  least  at  solstice.  Results  for  the  equinox  are  more  equivocal. 
Multiple  Layer  Model 

All  of  the  infrared  calculations  made  thus  far  have  been  bared  on  a 
two-stream  approximation  with  absorption,  emission  and  scattering.  Analytic 
solutions  within  layers  with  constant  ratios  of  scattering  are  matched  at 
the  boundaries  to  give  the  constants  in  the  equations  for  the  net  and  total 
fluxes.  The  equations  for  these  constants  are  of  the  form 


Ki  '  “iJi+si 


Ki  ■  Vi-i+BiKi.  1*1 

Ji  ‘  DiJi-i+EiKi-i+Fi 
JN  ' 


2  through  N 


An  algebraic  solution  was  made  for  N  =  3,  but  the  number  and  complexity  of 
the  terms  made  it  apparent  that  direct  solutions  for  larger  N's  were  imprac¬ 
tical,  while  matrix  inversion  techniques  were  time-consuming  and  involved 
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serious  Problem  with  small  differences  between  large  terms.  The  approach 
^  CUrre"tly  "°rked  «  involves  aerri„8  „p  the  e,uatIo„e 
Ki  "  JiYOO+e(i) 

setting  up  recursion  formulae  for  Y,  „  s,  and  ,,  a„d  solvlng  £or  ^ 

thC  C’Uatl0n  *”  ^  ~  *  «-  coefficients  allows  cancella- 

tion  of  many  large  terms,  and  the  final 

final  equations  appear  as  if  they  should  be 

computationally  stable.  Pinal  Caching  of  these  equations  and  programing  are 

““  C°mPleted  -11  -p  f  5  separate  cloud  layers  with 

UP  to  3  different  sets  of  optical  parameters. 
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APPENDIX  I 


THEORY,  INSTRUMENTATION,  CALIBRATION  AND  EXPERIMENTAL 
PROCEDURE  FOR  DETERMINING  OPTICAL  EXTINCTION  COEFFICIENTS 
ARISING  FROM  ABSORPTION  AND  SCATTERING  FROM  ATMOSPHERIC  AEROSOLS 

G .  Shaw 

DESCRIPTION  OF  EXPERIMENTAL  METHOD 

The  direct  solar  radiation  at  wavelength  X,  I0(X),  decreases  in  mag¬ 
nitude  during  its  progressive  entry  into  the  atmosphere  because  of 
scattering  and  absorption  by  atmospheric  gases  and  particulates  (aero¬ 
sols).  In  the  visible  and  near  infrared  part  of  the  spectrum,  absorption 
by  atmospheric  gases  is  primarily  due  to  excitation  of  vibration  rotation 
bands  of  tri-atomic  molecules  such  as  carbon  dioxide,  water  vapor  and  ozone. 
Non-resonant  scattering,  on  the  other  hand,  occurs  from  all  molecules 
and  particulates.  Resonant  scattering  from  atoms  and  molecules  and 
emissive  processes  such  as  Raman  or  fluorescent  scattering  can  be  neglected 
for  the  experiment  reported  here. 

We  make  the  assumption  that  the  Lambert-Beer  law  of  atmospheric 
attenuation  is  valid,  and  therefore  the  depletion  in  intensity  arising 
from  absorption  and  scattering  in  a  height  increment  dh  is  expressible 
as, 

dl(A,h)  =  -  I(A,h)  sec  (z)  dh  B(A,h)  (1) 

where  I^A,h)  ■  direct  solar  intensity  at  wavelength  A  and  at  geometric 
height  k 

z(h)  -  zenith  angle  of  incoming  solar  ray  at  height  h 
B(A,h)  -  volume  extinction  coefficient  (km  ^). 

Integration  of  eq.  1  leads  directly  to  an  expression  for  the  solar  in¬ 
tensity  at  h  lght  h  as, 

58 

A1 


(2) 


I(A,h)  -  I0 (A)  EXP 


EXP|-sec(z)  J 


8(h)dh 


MA>  '  8°Ur  lnten”lty  top  layers  o,  the  at™- 


sphere. 


The  optical  depth  T<X,h)  l8  defined  ae 


T(l. h)  -  f  6(h)dh 


It  le  assumed  that  a  t  a(h,  which  la  equlvale„t  to  neglecting  the 

tore  of  the  earth  or,  said  another  way,  the  atmosphere  Is  approxi¬ 
mated  as  being  place  parallel,  this  assumption  Is  valid  only  provlded 
that  the  solar  elevation  angle  BO*-,.  Is  target  than  about  »\  For 
lower  elevation  angles  the  Bemporad  function  Can  be  used. 


The  approximation  to  the  equation  of  transfer  represented  by  eq  1 
can  become  Invalid  for  atmospheres  containing  an  emission  term  (such 
ss  fluorescent  scattering,  or  Soever  multiple  scattering  or  large 
optical  depths  are  encountered.  For  radiation  In  the  visible  and  near 
Infrared  region  passing  through  the  earths  atmosphere,  however.  e„.  1 

represents  a  sufficiently  good  degree  of  approximation  for  clear  atmo¬ 
spheres  (with  ,W«1>.  The  requirements.  (hen>  for  eqs  q  ^ 

1  ■  .  . 


»S  valid  are:  (1)  the  atmospheric  optical  depth  must  be  small  (,<„, 
(2)  the  solar  elevation  angle  must  be  greater  than  about  20*.  I„ 


addition,  equations  are  Invalid  for  wavelength  regions  where  there  Is 
strong  absorption  by  gaseous  components  or  resonant  scattering. 


The  expression  provided  by  ea  1  f nrBC  , 

oy  eq.  i  forms  the  basis  for  the  inter- 


pre cation  of  the  optical  airborne  equipment  that  Is  the  subject  o,  this 


report.  Direct  solar  inte  isif-v  at-  „  , 

y  8everal  narrow  wavelength  intervals 


AA  - wavexengen  intervals 

X  '  '°2>  ""  "°"It0r0d  “1C"  “  airborne  photometer.  Changes  In  al- 

f  4  II 


tltude,  dh.  resulted  1„  small  changes  In  the  solar  Intensity  1„  accor¬ 
dance  with  the  expression  given  by  eo  1  TC  l 

ven  b)  eq.  1.  Thus,  by  noting  the  magnitude 


of  the  solar  Intensity,  I. (A)  and  the  differential  change  In  intensity 
vith  height,  dI(A ,h)/dh,  it  was  possible  to  quantitatively  evaluate 
the  optical  extinction  coefficient,  8(A,h),  by  using  eq.  1.  Furthermore, 
by  employing  multiple  wavelengths,  one  can  then  evaluate  the  ex¬ 
tinction  coefficient  at  multiple  wavelengths.  Knowledge  of  6(l.h)  at 
several  diverse  wavelengths  enables  one  to  deduce  some  parameters  that 
describe  the  site  spectrum  of  the  atmospheric  aerosols  contained  in  the 

neighborhood  of  h. 

The  evaluation  of  the  optical  parameter  »(X.b>  is  of  fundamental 
interest  for  modeling  the  radiative  transfer  processes  in  the  atmosphere, 
in  particular,  our  surveillance  of  the  size  distribution  and  height  dis¬ 
tribution  of  the  atmospheric  aerosols  has  important  significance  with 
regard  to  assessing  the  "radiation  climate"  over  the  pack  ice  and  over 

the  tundra  plains  of  northern  Alaska. 

THE  AIRBORNE  PHOTOMETER- INSTRUMENTAL  DETAILS 
An  airborne  photometer  was  constructed  for  the  optical  study  and  is 
illustrated  in  Figure  1.  The  photometer  is  hand-held  and  is  kept  pointed 
directly  at  the  sun  by  use  of  a  peep  sight.  The  photometer  is  csrried  sloft 
normally  by  a  Cessna  180  aircraft!  the  solar  disk  was  viewed  by  pointing 
the  instrument  through  a  five  centimeter  hole  in  the  aircraft  window. 

An  interference  filter  is  used  to  define  a  narrow  wavelength  region 
(AX/X  -  0.02) ;  provisions  are  made  to  change  the  filter  so  that  several 
wavelengths  can  be  used.  A  FIN  doped  silicon  photodetector  wit'  an 
active  area  of  1  cm2  is  employed  in  the  photometer;  it  has  a  usable 
response  extending  from  about  0.38  microns  to  1.15  microns.  Ihe  <»- 
coming  light  beam  is  modulated  with  a  mechanical  chopper  wheel  that 
allows  the  highly  temperature-dependent  diode  leakage  current  to  be 
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Gain  Selection 


to  Chart  Recorder 


Figure  1.  Illustration  of  portable  airborne  sun-photometer. 
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separated  from  the  component  of  electric  current  that  arises  from  light 
impinging  upon  the  photodetectoi . 

The  electronic  circuitry  is  designed  to  amplify  the  alternating 
electric  current  from  the  photodiode  to  a  level  sufficient  to  operate 
a  portable  chart  recorder.  Operational  amplifiers  with  controlled 
feedback  are  used  in  the  electronics  to  insure  constant  amplification 
under  a  wide  range  of  environmental  conditions.  The  schematic  diagram 
of  the  photometer  electronics  is  shown  in  Figure  2. 

The  solar  intensity  as  observed  near  the  ground  undergoes  a  con¬ 
siderable  variation  in  magnitude  as  solar  elevation  angle  and/or  wave¬ 
length  are  changed,  and  therefore  it  is  mandatory  to  incorporate  several 
amplification  levels.  Four  discreet  gain  levels  are  available,  eich 
^^er^n8  by  a  factor  of  five,  and  these  enable  one  to  obtain  a  usable 
signal  over  a  dynamic  range  of  approximately  10^. 

The  chart  recorder  and  photometer  both  operate  from  a  battery  pack; 
the  entire  system  weighs  approximately  10  kg  and  occupies  about  0.25  m3. 

FLIGHT  DESCRIPTIONS 

As  mentioned  previously,  the  direct  sun  photometer  was  flown  aboard 
an  aircraft  (Cessna  180).  All  flights  were  made  near  the  vicinity  of 
Barrow,  Alaska,  and  started  at  the  Naval  Arctic  Research  Laboratory.  A 
total  of  eight  flights  were  conducted  from  the  period  4  April  1972  to 
18  April  1972.  The  first  half  of  April  was  chosen  because  of  favorable 
chances  of  clear  ski.es  and  also  because  concurrently,  surface  studies 
were  carried  out  further  north  on  the  pack  ice  at  the  A'.DJEX  camp. 

Table  1  presents  a  brief  synopsis  of  all  eight  flights: 


62 

A4 


CHOePER 


COROER 


TABLE  1 


FLIGHT  SCHEDULE  -  BARROW,  ALASKA 


Flight  No. 

Date 

Time  (GMT) 

Max  Altitude 

(1000  ft) 

1 

4-  6-72 

0151-0251 

10 

2 

4-  6-72 

2120-2400 

10 

3 

4-  7-72 

1930-2014 

10 

4 

4-  8-72 

2025-2127 

10 

5 

4-11-72 

0038-0153 

12 

6 

4-12-72 

0118-0304 

15.5 

7 

4-16-72 

2201-2238 

10 

8 

4-18-72 

0037-0120 

10 

A  typical  flight  would  involve  climbing  at  a  constant 

rate  to  some 

specified 

altitude  ceiling,  normally  taken  as  10,000  feet,  all 

the  time  keeping  the 

instrument  directed  toward  the  sun  and  recording  the 

output  on  a  chart  recorder 

As  500  ft  increments  in  altitude  were  passed,  marku  were  placed  on 

the  chart 

record .  Temperature,  sky  conditions  and  other  relevant  information  were  logged 
during  the  flights. 

In  order  to  keep  the  instrument  pointed  at  the  sun  through  the  aper¬ 
ture  in  the  aircraft  window,  it  was  necessary  to  fly  at  a  constant  azimuth 
approximately  +40°  from  the  sun's  vertical.  This  was  because  of  the  geo¬ 
metry  imposed  by  the  overhanging  aircraft  wing  and  wing  strut  which  would 
occult  the  solar  disk  if  the  azimuth  was  not  properly  adjusted.  A  typical 
flight,  because  of  the  constraints  imposed  on  bearing  and  climb-rate,  would 
extend  out  some  100  miles  from  the  Naval  Arctic  Research  Laboratory  at 
Barrow.  Unfortunately,  the  assumption  of  horizontal  homogeneity  over  dis¬ 
tances  of  this  scale  is  not  necessarily  valid,  especially  during  certain 
synoptic  conditions  such  as  when  a  frontal  system  is  located  nearby. 

To  check  against  the  possibility  of  varying  airmasses,  we  would  occasionally 
maintain  a  constant  altitude  for  one  or  two  minutes  <>nd  check  to  insure 
that  the  solar  intensity  remained  constant.  Upon  occasion,  we  found,  by  using 
this  method,  that  the  solar  intensity  would  inexorably  decrease  (or  increase) 
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to  indicate  that  we  were  entering  into  a  more  (or  less)  turbid  air  mass. 

The  occurrence  of  thin  cirrus  clouds  aloft  was  of  continual  con¬ 
cern.  In  fact,  at  certain  times  when  the  sky  was  apparently  clear,  the 
solar  intensity  would  begin  to  decrease  while  flying  at  a  constant  altitude 
until,  eventually,  one  would  become  aware  of  very  tenuous  cirrus  ice 
clouds  in  the  vicinity  of  the  sun.  In  some  instances,  the  thin  cirrus 
striations  or  other  horizontal  homogeneities  invalidated  the  results 
of  a  flight.  Also,  instrumental  problems  occasionally  occurred,  notably 
freezing  of  the  chopper  motor,  with  an  attendent  noisy  signal. 

An  inspection  of  the  analog  records  and  log  book  indicates  that 
usable  information  can  be  obtained  from  four  flights.  The  data  accu¬ 
mulated  on  other  flights  is,  for  one  reason  or  another,  suspect;  therefore, 
these  records  will  not  be  analyzed. 

During  a  typical  flight,  the  optical  filters  were  placed  in  the 
photometer  one  at  a  time  and  data  would  be  accumulated  for  each  color 
by  ascending  or  descending  through  the  absorbing  and  scattering  layers. 
Because  of  the  relatively  long  time  involved  in  reaching  a  ceiling  ele¬ 
vation  of  10,000  ft  (this  would  take  about  a  half  hour)  we  limited  the 

O 

number  of  filters  used  to  two  colors  centered  at  5000  and  7000  Angstroms 
although  a  small  amount  of  information  was  also  accumulated  at  a  wave- 

O 

length  of  4000  Angstrons. 

DATA  REDUCTION 

The  data  consists  of  analog  chart  recordings  of  photometer  voltage 
as  a  function  of  altitude,  and  also  miscellaneous  notes  on  meteorological 
parameters  such  as  sky  appearance  and  temperature.  Supplementary  infor¬ 
mation  consi  ts  of  synoptic  weather  maps  at  the  surface  level  as  well  as 
at  850  mb,  700  mb  and  500  mb  levels. 
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Calibration  Procedures 


In  order  to  interpret  the  chart  recordings  one  must  calibrate  the 
photometer.  Calibration  of  the  photometer  is  necessary  for  several 
reasons.  Firstly,  one  must  insure  that  the  relationship  between  the 
output  signal  and  the  incoming  light  intensity  is  linear.  Secondly, 
the  relation  between  the  output  voltage  and  input  intensity  in  absolute 
units  is  desirable  because  it  allows  one  to  assess  some  additional  in¬ 
formation  about  the  optical  depth  of  the  atmosphere. 

The  transfer  function  (ratio  of  output  voltage  to  input  electric 
current)  for  the  electronic  section  of  the  photometer  was  determined 
experimentally  by  replacing  the  photodiode  by  an  alternating  current 
source.  The  electronics  yielded  a  linear  response  over  a  considerable 
range  of  input  current.  The  linearity  of  the  photometer,  as  well  as  an 
estimation  for  the  absolute  calibration  constant,  was  further  obtained 
by  utilizing  a  500  watt  quartz  halogen  lamp  as  a  source  of  known  radiation. 
A  range  of  input  intensities  covering  several  decades  was  obtained  by 
using  calibrated  neutral  density  filters.  By  this  technique  it  was 
determined  that  the  photometer  is  linear  to  an  accuracy  of  better  than 
two  percent  over  a  range  of  3  decades  of  input  light  intensity. 

The  output  voltage  of  the  photometer  that  would  be  developed  for 
zero  air  mass  is  of  special  importance.  Knowledge  of  this  constant 
(for  each  filter  used)  enables  one  to  specify  immediately  the  total 
vertical  extinction  or  optical  depth  by  use  of  eq.  2. 

Successful  application  of  this  method  depends  upon  knowledge  of 

the  magnitude  of  the  solar  intensity,  I0(A),  incident  upon  the  atrao- 

_2 

sphere.  Fortunately,  the  solar  intensity  in  absolute  units  (pw  cm 
nm”l)  has  now  been  determined  to  an  estimated  accuracy  of  several 
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percent  over  the  entire  visible  region  (Thekaekara,  1968)  and  tabulated. 
Knowledge  of  the  solar  topside  intensity,  Io(X),  and  the  absolute  in¬ 
strumental  constant  k(X)  then  allows  one  to  a-priori  specify  the 
optical  depth  of  the  atmosphere  (viz  eq.  2)  by 


t (h,X) 


-1 

sec(z) 


£n( 


■V(X,h). 
V0(X)  ‘ 


where  V0(X)  =  k(X)I0(X) 


(4) 


k  =  absolute  instrumental  calibration  constant  (volts  watt 

2. 

nm  cm  ) 

V0  =  zero  air  mass  voltage. 

As  a  final  point,  it  is  noted  that  the  optical  depth  of  the  atmo¬ 
sphere  can  be  determined  in  an  independent  fashion  by  using  the  Langley 
method  of  curve  fitting.  In  order  to  illustrate  this,  we  note  that  a 
plot  of  log  V(XQ ,h)  against  sec(z)  would  by  eq.  4  be  linear.  Further¬ 
more,  the  slope  of  the  linear  relation  is  proportional  to  optical  depth, 
t(X),  and  the  intercept  of  the  line  onto  the  ordinate  will  yield  the 


numerical  value  for  V0(X)  corresponding  to  topside  incident  solar  radiation 


I0(X)  can  be  inferred  in  an  independent  fashion  from  the  method  utilizing 
direct  calibration  from  a  standard  lamp.  The  two  independent  methods  of 

O 

determining  VD  for  X=5000  A  agreed  to  within  6%:  this  is  considered  to 
be  within  the  limit  of  the  standard  lamp  irradiance  and  might  be  improved 
upon  in  a  m<  exhaustive  study. 

Inference  of  Total  Extinction  Coefficient 

The  extinction  coefficient,  6(h,X),  may  be  obtained  from  the  data 
by  use  of  eq.  1.  In  order  to  evaluate  solar  zenith  angle  at  the  time 
of  observation,  the  celestial  spherical  triangle  having  vertices  at 
the  solar  disk,  the  local  zenith,  and  at  the  celestial  pole,  is  solved 


to  yield 
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(5) 


cos(z)  ■  sin(£)sin(<5)  +  cos(i)cos(6)cos(T  ) 

o 

where  i  =  local  latitude 

6  ■  solar  declination  .at  time  of  observation 

Tq  =  hour  angle  referenced  to  solar  noon 
The  hour  angle  of  the  sun,  T,  is  determined  by, 

Tq  =  LST-SLD-LLD-E  (6) 

Where  LST  =  local  standard  time  in  radians 

SLD  =  standard  longitude  difference  =  10  hr 
LLD  ■  local  longitude  at  observation  place 
E  =■  equation  of  time 

For  the  wavelength  of  interest,  the  extinction  coefficient  g(h,A) 
is  composed  of  additive  terms  arising  from  absorption  by  atmospheric 
ozone,  3  (h,A),  Rayleigh  scattering  by  atmospheric  molecules,  3  (h,A) 
and  by  absorption  and  scattering  from  the  aerosols,  3D(h,A).  Written 
mathematically 

3(h,A)  “  3g(h,A)  +  3R(h,A)  +  3D(h,A)  (7) 

Rayleigh  Scattering 

The  extinction  term  arising  from  Rayleigh  scattering  is  propor¬ 
tional  to  the  number  density  of  molecules  and  hence  to  atmospheric 
pressure.  Furthermore,  the  Rayleigh  scattering  term  has  a  A-^  depen¬ 
dency  and  so  we  can  write 

6R(h.X)  -  6r («.».)  £$■  (i*-/  (8) 

where  P(h)  =  atmospheric  pressure  at  height  h. 

Taking  A „ “5000  A  6R(o,A„)“l. 716x10  km  A  (Elterman,  1968).  The  varia¬ 
tion  of  pressure  with  height  is  taken  from  the  Standard  Atmosphere  tabu¬ 
lations  (Valley,  1965). 

Atmospheric  Ozone 

Atmospheric  ozone  has  broad  absorption  features  in  the  mid-visible 
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rche  Chappius  bands)  and  hence  its  influence  on  our  data  must  be  con¬ 
sidered.  Figure  3a  illustrates  the  molecular  absorption  cross-section 
for  ozone  as  a  function  cf  wavelength.  Figure  3b  presents  the  mean  sea¬ 
sonal  and  latitudinal  variation  in  columnar  ozone  amount  in  atmo-cm. 

The  majority  of  atmospheric  ozone  is  contained  in  the  lower  stratosphere 
(Dutsch,  1971)  and  the  optical  depth  arising  from  overburden  of  ozone 
is  constant  for  the  altitude  range  of  the  Cessna  180  flights  at  Barrow; 
however,  some  correction  of  ozone  optical  depth  with  altitude  is  made 
for  the  radiometer  data  acquired  on  the  later  series  of  Convair  990  flights 
that  extended  to  38,000  feet.  For  these  cases,  the  vertical  distribution 
of  ozone  was  scaled  from  data  presented  by  Dutsch  (1971)  for  the  appro¬ 
priate  season  and  latitude.  It  perhaps  should  be  mentioned  that  the  verti¬ 
cal  overburden  of  atmospheric  ozone,  as  can  be  seen  from  inspecting  Figure 
1  (Section  V),  takes  on  its  largest  seasonal  value  in  spring  at  high  latitudes 

fold  thus  the  absorption  effects  from  ozone  for  our  Barrow  data  obtained  in 
April  are  quite  pronounced. 
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molecular  absorption  coefficient 


WAVELENGTH  (A) 


J 
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•  * 

Figure  3a.  Molecular  cross-section  of  ozone  as  a  function 
of  wavelength. 


MONTH 

Figure  3b.  Mean  seasonal  and  latitudina'  variations  in 
columnar  atmospheric  ozone  amount. 
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